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I. INTRODUCTION 



The A resonance plays a prominent role in intermediate energy nuclear physics [|I|. 
Experimentally, it can be seen as a clear resonance of the partial amplitude f2i2j{^)i 
the I212J = -fss channel in t^N scattering experiments, at a Center of Mass (CM) energy 
^/s = Ma = 1232MeV and with a width Fa = 120MeV 0. In this paper we are concerned 
with the possibility of generating this resonance using known information at threshold, chiral 
symmetry constraints and exact unitarization of the amplitude. 

The practical description of resonances requires, in general, the use of some unitariza- 
tion method since exactly at the resonance energy the amplitude becomes purely imaginary 
and takes the maximum value allowed by unitarity [Q. On the other hand, most unitariza- 
tion methods are based on perturbation theory and hence require a resummation of some 
perturbative truncated expansion of the full amplitude. This obviously requires the exis- 
tence of an energy region where the unitarized amplitude approximately coincides with the 
perturbatively expanded one. This region is typically close to threshold, since exactly at 
that point the amplitude is purely real, and unitarity sets no constraints on it. In other 
words, the amplitude is justifiedly unitarizable if perturbation theory works reasonably well 
somewhere, for instance at threshold. If this is not the case there is nothing much one can 
do about it; one can formally unitarize the amplitude but the corresponding predictions at 
threshold will be in general very different from those obtained in perturbation theory. 

The modern way to effectively incorporate chiral symmetry and departures from it, is 
by means of Chiral Perturbation Theory (ChPT) [Q. Being a perturbative Lagrangian ap- 
proach it preserves exact crossing, perturbative renormalization and unitarity and allows for 
a bookkeeping ordering of Chiral symmetry breaking. All detailed information on higher en- 
ergies or underlying microscopic detailed dynamics is effectively encoded in some low energy 
coefficients which, for the time being, can only be determined experimentally. For processes 
involving only pseudoscalar mesons the expansion parameter is taken to be p^/(47r/)^ with 
p the four momentum of the pseudoscalar meson and / the weak pion decay constant [Q. 
Thus, this expansion works better in the region around threshold. Away from threshold, 
however, the expansion breaks down and the violation of unitarity becomes increasingly 
strong. Recently, the use of unitarization methods for pseudoscalar mesons has been shown 
to work, i.e. the predicted unitarized amplitudes reproduce the threshold region, and provide 
definite theoretical central values and error estimates of the phase-shifts away from threshold 
and up to about 1 GeV [^]. 

The extension of ChPT to include both mesons and baryons as explicit degrees of freedom 
becomes possible if fermions are treated as heavy particles but in a covariant framework 
0, yielding to the so called Heavy-Baryon Chiral Perturbation Theory (HBChPT) Ij^-Q. 
Here, the expansion to order = 1, 2, 3, ... is written in terms of a string of terms of the 
form e^/(/2'M^+^-2'), with / = 1, . . . , [(A^ + l)/2]. Here, / and M are the weak pion 
decay constant and baryon mass respectively. The quantity e is a generic parameter with 
dimensions of energy built up in terms of the pseudoscalar momenta and the velocity f ^ 
(f ^ = 1 ) and off-shellness k of the baryons, being the latter defined as usual through the 

o o 

equation pb =M v + k, with ps the baryon four momentum and M the baryon mass up 
to corrections generated by higher orders in the HBChPT expansion. Practical calculations 
|T^-|r^ show, however, that the convergence rate of such an expansion may not be as good 



1 



as it was in the purely mesonic case. Even at threshold, where the finite baryon mass efi^ects 
should be minimal, there appear sizeable corrections to the scattering lengths for the lowest 
partial waves. The situation obviously gets worse as one departs from the threshold region. 
In the particular case of ttN scattering, a systematic calculation has been done, so far, up 
to third order |[Tl],|12[, i.e. up to and including terms of order 1//^ (first order), l/(/^M) 
(second order) and 1//^ and l/(/^M^) (third order). A fit to threshold properties allows to 
extract the low energy constants, with the result that the second order contribution turns out 
to be larger than the first order one!0. For instance, in the P33 channel Mojzis obtains 
that the scattering length Ogg changes from 34.5 GeV~^ at first order, to 80.5 GeV~^ at 
second order and to 81.4 GeV~^ at third order. Obviously, any attempt to unitarize this 
amplitude based on the smallness of the second order term with respect to the first order 
one, will formally reproduce HBChPT but it will likely fail numerically to reproduce the 
corresponding phase shifts even in the vicinity of the threshold region. Thus, to be consistent 
with the calculation of Ref. ITTfl the used unitarization method should take into account the 



slow convergence rate of HBChPT, and ensure, to start with, that it numerically reproduces 
the scattering data around threshold. 



The A resonance can be included as an explicit degree of freedom within HBChPT |13 - 15 



This requires the introduction of new parameters into the Chiral Lagrangian. On the other 
hand, the unitarization of the amplitude via the Inverse Amplitude Method (lAM) does not 
introduce new parameters and has recently been proposed [l^, but the second order terms 
have been considered to be small. As a consequence, the low energy constants turn out to 



be very different from those found within HBChPT [|ll|] (see also the discussion below). In 



particular, their value 619 = —23.13, implies a pion nucleon coupling constant Qt^nn = 19.23, 
a rather odd one, compared with the experimental one QirNN = 13.4 ± 0.1. In the present 
work we show how the lAM should be modified in order to simultaneously i) implement 
exact unitarity , ii) comply with HBChPT at threshold and iii) describe the A resonance 
without having to introduce additional explicit parameters than those needed at threshold. 



II. PARTIAL WAVE AMPLITUDES 



We rely heavily on the notation of Ref. [|TT| and refer to that work for more details. The 
TcN scattering amplitude is given by 



T^N = HM V + p', a') A^ + B 



(1) 



where (M v + p,q) and (M v + p', q') are the incoming and outgoing nucleon and pion 

O 

momenta respectively. The bare nucleon mass M and the velocity f ^, ( f ^ = 1 ) determine 
the baryon off-shellness p and p' and a and a' are the spin indices. The supperscrit "±" 



^Both Ref. |jTl|] and Ref. |jT^ make a fit of the low energy contribution either by extrapolation 
of close to threshold data to threshold or by extrapolation of the theory to the close to threshold 
region respectively, but only the first one displays the contributions of the several orders explicitly, 
so we will mainly refer to Ref. in this regard. 
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corresponds to the isospin decomposition T'"^ = T"'"(5^" + ie''°^T~rc, with a and h the incoming 
and outgoing pion isospin states in the cartesian basis respectively. and 5^ are scalars 
depending on the Mandelstam variables s, t and u which for on-shell pions and nucleons 
verify s + t + u = 2{m? + M^), with m the physical pion mass and M the physical nucleon 
mass. In the CM frame, the partial wave amplitudes are written as (the subscript ± stands 
forj = /±l/2), 

fUs) = e^^^lW^^^ = ^ ( (E + M)[AHs) + (v^ - M)BHs) 

-{E- M)[Af^,{s) -{^s + M)S±i(s)]) (2) 



where, for elastic scattering, q = \q\ is the CM momentum and E = y'g^ + M'^ and u = 
+ the CM nucleon and pion energies respectively. Finally, y/s = E + u is the total 
CM energy and S^{s) are the corresponding phase-shifts. The projected amplitudes Ai{s) 
and Bi{s) are defined by the integrals 

{aHs),BHs)) = j'^dz{A^{s,z),B^{s,z))Pi{z) (3) 
with z = cos6 and 6 the CM scattering angle. In addition, A^ and B^ can be written as 

A^ = {a^ + (4) 

5" = (-^ + 4V^ (5) 



4M 

with t = — 2g^(l — cos6) and u = + — 2Euj — 2q^cos9. An expansion of both and 
along the lines of HBChPT has been given in Ref. [^] up to third order, i.e. including 
terms of order e/ p (first order), e^/{pM) (second order) and e^//^ and e^/(/^M^) (third 
order). There, however, in some cases the full nucleon mass dependence has been retained. 
This makes the discussion on orders a bit obscure^. We have preferred to further expand 
any observable in terms of 1/M or 1//^. As a consequence and in contrast to Ref. [0, there 
is no dependence of D,F, and higher partial waves on the low energy constants. Thus, the 
low energy parameters are solely determined by S and P partial waves. One thus gets the 
following expansion, 

/,4=/,'r+/,'r +/,?*+■■• (6) 



where 



/,'r = ^4'/>*0 (7) 

A2)± _ m (l,2)±/^\ /„v 



A3) ± _m (3,3) ± / w \ m (1,3) ± / w N 



^For instance, the P33 scattering length contains some corrections m/M to all orders ( see eqs.(90- 
92) in that reference ), instead of neglecting all pieces of order l/(/^M^), l/(/^M) or higher as 
dictated by the expansion assumed in the Chiral Lagrangian. 
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^(ra,m)± dimensionless functions of the dimensionless variable cu/m, independent of /, M 
and m, whose analytical expressions are too long to be displayed here. The convenience of 
the double superscript notation will be explained below; n + 1 indicates the power in 1// and 
m the total order in the HBChPT counting. The unitarity condition Im/;"^^ = —q becomes 
in perturbation theory 

ImtS^/)± = Imti'f^ = Imtp/)^ = (10) 
Imf/)^ = ^|4^/)^r (11) 

The scattering lengths af^^ and effective ranges bf^^ are defined by 

Re/±± = g2^(a±± + g26±^ + ---) (12) 

Obviously, the expansion of Eq.(^) for the amplitudes carries over to the threshold parame- 
ters, af^ and (see below). 



III. RE-FITTING THE LOW ENERGY CONSTANTS 



As we have said, in Ref. |]TT|, a systematic expansion for the fixed-t amplitudes and (3^ 
was undertaken, but some higher order terms were retained in the partial wave amplitudes. 
This makes the discussion about orders a bit unclear, since different orders are mixed. In 
all our following considerations we further expand the amplitudes, and consequently the 
threshold and close to threshold parameters in the spirit of HBChPT. For the scattering 
lengths we have 

a',,.,nr-» = ^„<M. + ^„<«. + + ^„,3,3,+ (13) 

and a similar expression for 62/ 2 j"^^'^^- With this expansion we may reanalyze the fit of 



Ref. ||TT|. For a better comparison we do so for the same experimental data, and with 
the same input values of the parameters, M = 939 MeV, m = 138 MeV, / = 93MeV and 
qa = 1.26. The result of the fit is presented in Table 1. As we see from the table the changes 
in the parameters are small, and are compatible within two standard deviations. This was 
to be expected since the difference between Mojzis's parameters and ours is of higher order. 
It is also noteworthy that the resulting = 10.68 is mainly made out of b^^ and a which 
contribute with 3.3 and 7.3 respectively to the total x^- The bad result concerning the 
cr-term is substantiated by the findings of Refs. HIDHT^ . 

After performing this re-fitting procedure we have found instructive to separate the 
contributions to the scattering lengths and effective ranges of the lowest partial waves, as 
done in Table 2. A clear distinctive pattern emerging from Table 2 is that the contribution 
of order l/(/^M^) is always rather small. Only in some cases, however, is the contribution 
of order 1//^ also small. This is so in the P33 channel in particular. Incidentally, let us note 
that for this channel the smallness of the third order contribution is due to the smallness of 
both the and the l/{pM) terms and does not stem from a cancelation between large 
contributions. So, it seems that close to threshold the 1/ P expansion converges faster than 
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the 1/M expansion. Actually, the terms in 1//^ and 1 / {pM) are of comparable importance. 
From here, it is clear that any unitarization method will only be consistent with the HBChPT 
approach, if it treats both the first and the second order as equally important. 





Mojzis (Ref. Q) 


This Work 




Set I 


Set II 


0.1 


fin 4- n IT? 

— Z.DU ± U.Uo 


— z. /z ± U.Uo 


02 


1.40 ± 0.05 


1.39+0.05 


OS 


-1.00 ± 0.06 


-1.11+0.06 


05 


3.30+0.05 


3.3+0.07 


61+62 


2.40+0.3 


1.8+0.4 


63 


-2.8+0.6 


-3.4+0.7 


66^ 


1.4+0.3 


1.9+0.3 


616 - 615 


6.1+0.6 


6.1 +0.6 


619 


-2.4+0.4 


-2.4+0.4 



TABLE I. Low energy constants obtained from fitting the threshold parameters a^_|_, &^_|_, cli± found 
in the present work and in Ref. p4[ |, the Goldberger-Treiman discrepancy and the nucleon cr— term to the 
HBChPT predictions. We take M = 939 MeV, m = 138 MeV, / = 93MeV and gA = 1.26. Experimental 
data are taken from Ref. |17|. The first column refers to the work of Mojzis jl^ , where the fixed t amplitudes 
and have been expanded using HBChPT but not the threshold parameters. In our case, we have 
further expanded the threshold parameters (see text for details) and fitted to the same data. The dependence 
of the Goldberger-Treiman discrepancy and the nucleon u-term on the fitted parameters can be found in 
Ref. |ll| . There is only one degree of freedom. Errors in any parameter have been calculated by changing 
the corresponding parameter until that the minimized with respect to the remaining parameters varies 
by one unit. The fit yields a x^i„=10.68 which is basically made out of and a which give 3.3 and 7.3 
contribution to the total respectively. The value of the total in Ref. Q is not given. 







l/(/'M) 


1/(/2m2) 


1//' 


Total 


Exp. (Ref. 0) 




-0.63 


-0.08 


+0.07 


-0.11 


-0.75 


-0.73 + 0.01 


o?,i 


1.27 


-0.36 


-0.07 


0.41 


1.25 


1.25 + 0.02 


"■3 3 


35.3 


48.1 


-1.6 


-0.3 


81.5 


81.5+0.8 


Ol3 


-17.6 


14.3 


-3.0 


-5.2 


-11.5 


-11.5 + 1.3 


O31 


-17.6 


11.8 


-1.7 


-9.6 


-17.1 


-17.2 + 0.6 


Oil 


-70.6 


86.9 


-3.6 


-43.1 


-30.4 


-30.4 +0.9 



TABLE II. HBChPT, lowest partial S and P wave- scattering lengths, a^j^j (in GeV"^'"^ units) for txN 
scattering in the in HBChPT, decomposed as a sum of terms of first order 1//^, second order, l/(/^M), and 
third order, l/{p]\,P) and !/(/"'). The sum of all the terms yields the total scattering length parameter. 
The experimental value obtained in Ref. is also quoted. 
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IV. UNITARIZATION METHOD FOR THE P33 CHANNEL 

Our unitarization method is based on assuming, as suggested by the perturbative calcula- 
tion, that the chiral expansion in terms 1//^, has a stronger convergence rate than the finite 
nucleon mass 1/M corrections. Such a hypothesis can only be supported by confrontation 



with experimental data, and indeed there are some recent theoretical attempts |18,19| to 



define a relativistic power counting not requiring the heavy baryon idea somehow retaking 



the spirit of older relativistic studies |]20|. Indeed, we show that a slightly modified version 
of the well-known lAM approach turns out to describe the P33 phase-shift satisfactorily, just 
using the low energy parameters determined at threshold, i.e., without re- fitting them in the 
unitarized case. The idea is to consider the expansion of the inverse amplitude in terms of 
m? / p. The expansion of the amplitude can be written as 



3 

fl^iuo, m, /, M) = j/^\uj/m , m/M) + !^t^^\uj/m , m/M) + . . . (14) 

where we have explicitly used that the functions ^(^"+1) are dimensionless, and hence depend 
only on dimensionless variables, such as u/m and m/M. For the sake of a lighter notation, 
the quantum numbers / / and J have been purposely suppressed. Perturbative unitarity in 
this expansion requires at lowest order, 

Im t(^) {u /m, m/M) =0 (15) 

lmt^^\uj/m,m/M) = — |tW(cu/m, m/M)P (16) 

m 

which in turn imply an infinite number of conditions in the 1/M expansion. The functions 
t^^^ and t'^^^ are only known in a further m/M expansion, 

t(2"+l)(^/rn,m/M) = t(2-+l,2n+l)^^/^) + ^t{2n+l,2n+2)(^/^) ^ ( ^(2„H-l,2n+3) 

(17) 

yielding Eq.(|^) and Eq.(|^) after a suitable isospin projection]^. For the inverse amplitude we 
get then 

1 p 1 t^^Huj/m, m/M) , , 



p3{s,m,f,M) mt^^\uj/m,m/M) [t^^\uj / m , m / M)Y 

Obviously, expanding in m/ p may be justified provided these corrections are small some- 
where. As we have shown, for the P33 channel, they are small precisely at threshold. So in 
the threshold region the unitarization scheme will, approximately, reproduce the perturba- 
tive result. At the same time, unitarity is exactly implemented. However, the 1/M terms are 
not small corrections at threshold, so we refrain from further "expanding the denominator" . 
In this way we keep the first mass corrections as equally important. The state of the art 



The relation between /^^_2J and /,^^ is given by j\^2i±\ = ft± ~ fi,2i±i = ft± + '^^i 
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in HBChPT calculations is such that only t^^'^\ t^^''^\ t^^'^^ and t^^'^^ are known. To ensure 
exact elastic unitarity we should keep only t^^'^^ in the denominator of the second term in 
Eq. (|18D above. This induces a tiny error thanks to the smallness of t*-^'^-* (correction 1//'*) 
in the P33 channel. Of course, it would be highly desirable to compute at least t^^''^^ and 
^(3,5) (^orders l/(/^M) and l/(/^M^) respectively) in order to be able to keep also t^^''^'^ in 
the denominator of the second term of the mentioned equation. After these remarks, we 
have 

1 P 1 



/3l3(s,m,/,M)|u„itarized m t^^ (u / m) + § t^^ {u / m) + (^)2 t{i,3)(^/^) 

t(3,3)(^/^) 

~'"[t(i'i)(a;/m)]2 ^^^^ 
At threshold, our formula yields a modified scattering length 

1 p 1 a(3,3) 

"^ r.,n.i^i2 (20) 



aunitarized m a^'D + ^ + (^)2 tt(l,3) [«(!,!)] 



tt(3,3) 



"^TTTTlm (21) 



aHPChPT - {^) a(^'^) 

which on view of Table 2 yields, for the P33 channel, the value agglunitarized = 80.2GeV~^, 
to be compared to a33|HBChPT = 81.0GeV~^, a compatible value with the experimental one. 
Notice that if we had expanded the denominator considering the second order contribution ( 
l/(/^M) ) to be small we would have obtained strictly the lAM method as used in Ref. ||16|| , 




aiAM m [a(i'i) 'M'[a(^'^)]^ ^M' [a^^^^)f ^M' [a^^^^)f ] [a^^^^)f 
yielding a3 3|iAM = 22.8GeV~^, a completely odd result. This explains why the low energy 



parameters recently found in Ref. ||I6| are so different from those found by Mojzis 
Actually, their value of 619 = —23.13 would yield a pion-nucleon coupling constant Qt^nn - 
19.2, completely out of question. 



V. NUMERICAL RESULTS 

From our previous discussion it is clear that at threshold our unitarized amplitude will 
reproduce very accurately and within error bars the HBChPT results and hence the exper- 
imental data. It is thus tempting to extend the P33 phase shift up to the resonance region 
and propagate the errors of the low energy parameters. We show in Fig.l and Fig. 2 the cor- 
responding phase shifts as the outcome of our unitarization method, Eq.( |1^), and compare 
them to the experimental vriV data 0. We use both parameter sets given in Table 1. In 
Fig. 2 we use the parameters determined in Ref. []rT| ( set I ) and in Fig.l those determined 
in the present work ( set II ). 
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1200 



1500 



1300 1400 
(MeV) 

FIG. 1. P33 phase shifts as a function of the total CM energy y/s. Experimental data are from Ref. 
The upper solid curve is the result of our approach (Eq. |l9|), with the corresponding propagated error bars 
(dashed-dotted curve). The intermediate solid curve represents the static limit prediction, M 00, also 
with error bands (dashed-dotted). The bottom solid line is the prediction of the conventional I AM method 
of Ref. In all cases, the low energy parameters are those labeled as set II in Table 1. The position and 
width of the resonance for the upper curve are Ma — 1267l25MeV and Fa = 213l54MeV respectively. 



As one can see from both figures, the prediction of Unitarized Heavy Baryon Chiral 
Perturbation Theory agrees rather well within errors with the data. It is fair to say, however. 



that for the parameters determined in Ref. |TT] the agreement seems much better, describing 
the data for larger CM energy values. We do not ascribe any particular significance to a 
very accurate agreement at larger energies since there is certainly more physics to be taken 
into account, like inelasticities. In both Fig.l and Fig.2, the error bars reflect only the 
uncertainty in the low energy parameters and have been obtained by means of a Monte 
Carlo gaussian sampling of the low energy parameters for any given CM energy value. The 
location and width of the A resonance come out to 



Ma 



1238: 



h22 



150+^? 



18 MeV (set I), 
MeV (set I), 



1267+1 MeV (set II) , (exp. 1232 ± 2 MeV) 
213tlt MeV (set II), (exp. 120 ± 10 MeV) 



(23) 



in agreement within two standard deviations with the experimental numbers, although with 
much larger errors. Given the good quality of our description from threshold up to the 
A resonance, one might even re-fit the parameters of our unitarized amplitude, and hence 
reduce the errors of the low energy constants entering the definition of the amplitude in the 
P33 channel. In this way we would reduce errors quoted in Eq. ( ^3]) as suggested in Ref. |]12 



^ We obtain Ma from the condition ^3 3 



{Ml 



tt/2, and the width Fa from 1/Fa 



Ma 



ds 



-Ml 
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1200 



1500 



1300 1400 
v/i (MeV) 

FIG. 2. The same as in Fig.l but with the low energy parameters labeled as set I. The position and 
width of the resonance for the upper curve are Ma = 1238l^gMeV and Fa — 150^3^ MeV respectively. 

To have an idea of the convergence rate of our calculation, we have also depicted in 
both figures the prediction in the static limit (M oo). As we see, there also appears 
a resonant behaviour, but at higher energies \fs ~ 1450MeV. Thus, though the bulk of 
the dynamics is contained in the static limit, the finite mass corrections, particularly the 
l/(/^M) contribution, are important to achieve an accurate description^. 

Finally we also show the results obtained within the conventional JAM approach, see 
for details Ref. |jT6|, for both sets of parameters of Table 1. As we already anticipated, the 
description is extremely poor mainly due to improper treatment of the l/(/^M) corrections. 
Thus within this framework, a fit to the data becomes possible only if huge changes in 
the low energy parameters are allowed. For instance, the authors of that reference get 
&16 ~ ^15 = 68.9, 6i9 = —23.13, 05 = 27.76, = 9.1, 03 = —8.7 ■ ■ ■ which strongly differ from 
the values of Table 1. 

The situation in other I212J channels might be somehow different and definitely deserves 
further study, though the results of the present paper regarding the P33 channel, will remain 
unchanged. We anticipate, as can be already seen from Table 2, that the agreement of 
our unitarizated threshold parameters with those stemming from HBChPT, would get a 
bit spoiled in some of the remaining channels, since the 1//^ is not a small contribution 
due to a, perhaps accidental, but effective cancelation of the leading, l/p, and the next 
to leading, l/(/^M), order terms. The situation might improve, however, if the low energy 
constants would be allowed to vary as to provide a better description of the scattering data 
from threshold up to the A resonance region. A systematic and detailed investigation of 
this topic will be presented elsewhere [pT|. 



^The l/(/^M^) correction turns out to be quite small. That was expected, since it neither provides 
a sizeable contribution at threshold (unlike the l/(/^Af) correction ) nor it is responsible for the 
restoration of unitarity (like the 1//^ correction) . 
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VI. CONCLUSIONS AND OUTLOOK 



We summarize our results. Heavy Baryon Chiral Perturbation Theory provides definite 
predictions for tfie ttN scattering amplitudes in the threshold region, but it violates exact 
unitarity if the perturbative expansion is truncated to some finite order. Hence it is unable 
to describe the A resonance in the P33 channel. The analysis up to third order shows 
that the leading finite nucleon mass correction, which is second order, is of comparable size 
to the static approximation and in fact it dominates the corrections at threshold. This 
sets conditions on the unitarization method suggesting an expansion in inverse powers of 
the weak pion decay constant but without making the heavy baryon expansion. Such an 
idea is supported by recent theoretical attempts to redefine a relativistic chiral counting for 
baryons. Our unitarization method is a slight but important variant of the lAM and provides 
a prediction for the phase shifts in the P33 channel, in terms of low energy parameters and 
their errors, as determined from HBChPT. The agreement with experimental data is good 
within uncertainties. It is also clear that with our unitarization method, the combination of 
parameters which enter the P33 amplitude, might be determined to a better accuracy than 
that obtained by only looking at threshold, since the central predicted values of the P33 
phase shift, fall mainly on top of the central experimental data. 
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